An easy method for finding many very large 
prime numbers* 



Leonhard Euler 



1. It suffices to explain this metliod with a single example of numbers con- 
tained in the formula 232aa + 1, where I shall investigate all the values of a in 
which this formula produces composite numbers; with these excluded, all the 
remaining numbers assumed for a will supply prime numbers. 

2. Now, whenever the formula 232aa + 1 yields a composite number, an 
equation like this will hold: 

232aa + 1 = 232xx + yy, 

because 232 = 8 • 29 is a numerus idoneus0 Therefore it will be 

232(aa — xx) ^ yy — 1. 

Now let us put y = 1 + 58z0 and dividing by 232 produces this equation: 
aa — XX — iz(29z + 1), or, because the unity needs to be taken both positive 
and negative, it will be 

aa — XX — — z(29z ± 1). 



•Presented to the St. Petersburg Academy on March 16, 1778. Originally published as 
Facillima methodus plurimos numeros primos praemagnos inveniendi, Nova acta academiae 
scientiarum Petropolitanac 14 (1805), 3-10. E718 in the Enestrom index. Translated from the 
Latin by Jordan Bell, Department of Mathematics, University of Toronto, Toronto, Ontario, 
Canada. Email: jordan.bell@gmail.com 

^Translator: A number n is a numerus idoneus (suitable number, convenient number, 
idoneal number) if it satisfies the following criterion. Let m be an odd number relatively prime 
to n such that for some relatively prime x,y, m = -\-ny^ . If the equation m = + "2/^ has 
only one solution with x,y >0, then m is a prime number. Thus if n is a numerus idoneus 
we can use it to find prime numbers m. Gauss: a positive integer n is a numerus idoneus 
if and only if for every quadratic form of discriminant —in, every genus consists of a single 
class, see Proposition 3.24 of David A. Cox, Primes of the form x^ +ny^ . It is also equivalent 
to h(— 4n) = 2''^^, where h is the class number and /i is defined in the following way. Let 
D = 0, 1 (mod 4) be negative and let r be the number of odd primes dividing D. If _D = 1 
(mod 4) then fi = r, and if D = —in, then /i is determined as follows: n = 3 (mod 4) then 
= r; n = 1, 2 (mod 4) then /i = r + l;n = 4 (mod 8) then ^ = r-|-l;n = (mod 8) then 
= r + 2. See Theorem 3.22 of Cox. For more about numeri idonei see Weil, Number theory: 
an approach through history, pp. 219-226. 

■^Translator: I don't see why Euler says y = 1 + 58z instead of j/ = 1 ± 582. y = 1 ± 582 
follows because 58 = 2 • 29, and both 2 and 29 are prime factors of the left hand side. 



1 



3. Now let us take in order all values 1,2,3 etc. for z and, because the 
formula ^^(292; ± 1) certainly has two factors and perhaps more, it can be 
represented as a product rs in one and perhaps more ways. Since it needs to be 
equal to the formula aa — xx = (a + x){a — x), let us put a + x = r and a — x = s, 
whence a = ^-j^, and the values of a to be excluded will be thus represented. 
Here it is evident that both factors r and s need to be both even or both odd. 

4. With this noted, let us take z = 1, and it will happen in two ways, because 
of the ambiguity of the sign, as either rs = 14 or rs = 15. Here the first value, 
since it is oddly even, is irrelevant to our purpose. On the other hand the latter 
value rs = 15 yields two exclusions, for it will either be r = 15 and s = 1 or 
r = 5 and s = 3, from which the values 8 and 4 need to be excluded for a. 

5. Now let 2; = 2 and it will be rs = 57 or rs = 59. The first admits two 
resolutions: 

r = 57, r = 19, 
s = 1, s = 3, 

the latter one: r = 59 and s = 1, from which these three exclusions arise: 
a = 29,a = ll,a = 30. 

6. Now let 2; = 3, and it will be either rs = 3 • 43 or rs = 3-44, whose 
resolutions happen like this: 



r = 129, r = 43, 
s = 1, s = 3, 



r = 66, r = 22, 
s = 2, s = 6, 



from which these four exclusions follow 



a = 65, a = 23, a = 34, a = 14. 

7. Let z = 4:, and it will be either rs = 230 or rs = 234. This numbers, 
since they are oddly even, give no exclusions. 

8. Now let z = 5, and it will be rs = 5 ■ 72 or rs = 5 ■ 73. The first of these 
yields the following resolutions: 

r = 180, r = 90, r = 60, r = 36, r = 30, r = 20, 
s = 2, s = 4, s = 6, s = 10, s = 12, s = 18. 

Then the values of a that need to be excluded are therefore 91, 47, 33, 23, 21, 19. 
The other value gives r = 73, s = 5 or r = 365, s = 1, from which the exclusions 
39 and 183 arise. 
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8aH Let z — 6, and it will be either rs = 3 ■ 173 or rs = 3 • 175, whence the 
following resolutions: 



r = 519, 173, 
s = 1, 3, 



r = 525, 175, 105, 75, 35, 25, 
s = l, 3, 5, 7, 15, 21. 



Therefore these exclusions arise: 



260, 



263, 



89, 



55, 



41, 



25, 



23. 



9. Now let z = 7, and it will be either rs = 7 ■ 101 or rs = 7 ■ 102. Here the 
first case alone gives these resolutions: 

r = 707, r = 101, 
s — 1, s — 7, 

and then the exclusions will be 354 and 54. So that we can prescribe some 
bounds for ourselves, we will omit in the sequel all exclusions greater than 300. 

10. Now let z = 8, and it wiU be either rs = 4 ■ 3 • 7 ■ 11 or rs = 4 • 233. Then 
the following resolutions arise: 



r = 462, 154, 66, 42, 
s = 2, 6, 14, 22, 



r ^ 466, 
s = 2, 



from which we conclude the following exclusions: 

a = 232, 80, 40, 32, 234. 

11. Let z — and it will be either rs = 9 • 130 or rs = 9 • 131. Then for the 
second 

r = 393, 131, 
s = 3, 9, 



one therefore excludes a = 198, 70. 

12. Let z — 10, and it will be either rs = 3 • 5 • 97 or rs = 5 • 17 • 17, then 



^Translator: The original has two sections 8. 
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r = 485, 291, 97, 
s = 3, 5, 15, 



r = 289, 85 
s = 5, 17, 



whence one excludes 

o = 244, 148, 56, 147, 51. 

13. Let z = 11, and it will be either rs = 3 • 11 • 53 or rs = 32 ■ 5 • 11, then 



r = 583, 159, 53, 
s = 3, 11, 33, 



440, 220, 176, 88, 80, 110, 44, 
4, 8, 10, 20, 22, 16, 40, 



hence 



a = 293, 85, 43, 222, 114, 93, 54, 51, 63, 42. 

14. Let z — 12, and it will be either rs = 6 ■ 347 or rs = 6 • 349, neither of 
which gives any exclusion. 

15. Let z = 13, and it will be either rs = 13 • 4 • 47 or rs = 27 • 13 • 7, then 



r = 94, 
s = 26, 



351, 273, 189, 117, 91, 63, 
7, 9, 13, 21, 27, 39; 



one therefore excludes 



60, 179, 141, 101, 69, 59, 51. 



16. The calculation can be continued on in this way as far as it seems fit. 
We shall represent the exclusions arising from all the values of z in the following 
table: 
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130, 141, 190, [202], 
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137, 143, 162, 169, 247, 271, 
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141, 145, 171, 287, 
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z 


Exclusions 


38 


212, 




39 


154, 


202, 


40 


156, 


174, 178, 242, 


41 


162, 


186, 


42 


160, 


220, 268, 


43 


164, 


[181], 195, 199, 211, 284, 


45 


208, 




46 


191, 


217, 257, 


47 


184, 




48 


185, 


241, 253, 


49 


[192], 260, 


50 


191, 


193, [215, 225, 279, 297], 


51 


196, 


236, 292, 


53 


202, 


223, 245, 260, 


54 


224, 




56 


216, 


240, 


58 


224, 


236, 296, 


59 


225, 


231, 233, 273, 


61 


239, 


241, 282, 


62 


237, 




63 


240, 




64 


246, 


282, 


65 


252, 




66 


274, 




67 


282, 




69 


263, 


265, 295, 263, 265, 


71 


278, 




73 


286, 




75 


286, 


290, 


78 


298. 





17. Among the numbers to be excluded for the letter a, some occur two or 
three times. This happens whenever the number 232aa + 1 has three or more 
factors; thus we see that the number 265 is excluded first by the value z — 18, 
then twice by the value z = 69. In fact in this case it turns out that the number 
232 • 265^ + 1 is made from the three factors 59 • 461 • 599, and even aside from 
this form the following three are recovered: 

1 ° 232 • 256^ + 10432, 

2 ° 232 • 342 + 40032, 
3° 232-352+40012, 

18. Let us also record in a table all the numbers that are to be excluded 
according to the following arrangement, and we note with an asterisk all those 
with occur more than once: 
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1 


11, 14, 19, 


2 


21 23* 25 [291 


3 


30, 32, 33, 34, 39, 


4 


40, 41, 42, 43, 47, 


5 


51*, 54*, 55, 56, [57], 58, 59, 


6 


60, 61, 63, 64, 65, 66, 68, 69, 


7 
8 


70, 75, 77, 78, 

80*, 83, 84, 85*, 88*, 89, 


9 


90, 91, 92, 93, 96, 97, 98, 


10 


101, 102, 103*, 105, 108, 109, 


11 


110, 111, 114, 116, [117], 


12 


120, 122*, 123, 125, 128, 129, 


13 


130, 131, 132, 133, 134, 135, 137, 139, 


14 


140, 141*, 143, 145*, 146, 147*, 148, 


15 


154*, 156, 


16 


160, 161, 162, 163, 164, [165], 166, 168, 169, 


17 


171, 178*, 179, 


18 


[181], 183, 184, 185, 186, 187, 


19 


190, 191*, [192], 193, 194, 195, 196, 198, 199, 


20 


202*, 203, 206, 207, 208, 209, 


21 


211, 212, 214, [215], 216, 217, 218, 


22 


220, 221, 222, 223, 224*, 225*, 


23 


231, 232, 233, 234, 235, 236*, 237, 239, 


24 


240*, 241*, 242, 243, 244, 245, 246, 247, 


25 


252, 253, 256, 257, 


26 


260*, 263*, 265*, 266, 268, 


27 


271, 272, [273], 274, 276, 278, [279], 
282*, 284, 286*, 287, 


28 


29 


290, 292, 293, 295, 296, [297], 298. 



Therefore with these numbers excluded, all the remaining written in place 
of a in the formula 232aa+ 1 produce prime numbers. The following table thus 
exhibits these values of a: 
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1, 


2, 


3, 


5, 


6, 


7, 


9, 


10, 


12, 


13, 


15, 


16, 


17, 


18, 


20, 


22, 


24, 


26, 


27, 


28, 


31, 


35, 


36, 


37, 


38, 


44, 


45, 


46, 


48, 


49, 


50, 


52, 


53, 


62, 


67, 


71, 


72, 


73, 


74, 


76, 


79, 


81, 


82, 


86, 


87, 


94, 


95, 


99, 


100, 


104, 


106, 


107, 


112, 


113, 


115, 



118, 119, 121, 124, 126, 127, 136, 138, 142, 144, 149, 

150, 151, 152, 153, 155, 157, 158, 159, 167, 170, 172, 

173, 175, 176, 177, 180, 182, 188, 189, 197, 200, 201, 

204, 205, 210, 213, 219, 226, 227, 228, 229, 230, 238, 

248, 249, 250, 251, 254, 255, 258, 259, 261, 262, 264, 

267, 269, 270, 275, 277, 280, 281, 283, 285, 288, 289, 

291, 294, 299; 

the last of these numbers comes to about 20 million0 



"Translator: If a = 299 then 232aa + 1 = 20,741,033. 
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